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Abstract 

We consider the problem of computing upper and lower bounds on the price of a 
European basket call option, given prices on other similar baskets. We focus here on 
an interpretation of this program as a generalized moment problem. Recent results by 
Berg & Maserick (1984), Putinar & Vasilescu (1999) and Lasserre (2001) on harmonic 
analysis on semigroups, the K-moment problem and its applications to optimization, 
allow us to derive tractable necessary and sufficient conditions for the absence of static 
arbitrage between basket straddles, hence between basket calls and puts. 

Keywords: Semidefinite Programming, Static Arbitrage, K-Moment Problem, 
Basket Options. 



1 Introduction 

We let p G R" +m , K G R n+m+1 , Wi G R n , % = 0, . . . , n + m and we consider the problem 
of computing upper and lower bounds on the price of an European basket call option with 
strike K and weight vector w : 

maximize/minimize po := E u (u)qX — Kq)+ , , 

subject to Ei„(wjx — Kj) + — p^ i = 1, . . . , n + m, 

with respect to all probability measures v on the asset price vector x G R/f , consistent with 
the (given) set of observed prices pi of options on other baskets. 

We implicitly assume that all the options have the same maturity, and that, without loss 
of generality, the risk-free interest rate is zero (we compare prices in the forward market). 
We seek non-parametric bounds, i.e., we do not assume any specific model for the underlying 
asset prices, our only assumption is the absence of a static arbitrage today (i.e. the absence 
of an arbitrage that only requires trading today and at maturity). 
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Here, we interpret (0) as a generalized moment problem. This approach was successfully 
used in [BP02J to get tractable bounds in dimension one and to show the NP-hardness of the 
multivariate problem (|TJ) . NP-hardness means that we have no chance of finding a direct and 
efficient method for detecting all arbitrage opportunities, here instead we look for a sequence 
of successively tighter price bounds. This means that outlandish arbitrage opportunities can 
be detected at little numerical cost while detecting finer price discrepancies has a higher 
theoretical complexity 

Recent results on multivariate moment problems (see |Sch91j . |PV99j or |CF00| ). semidef- 
inite programming (see |NN94j . |VB96j and |Nes00j ) and harmonic analysis on semigroups 
(see [BCR84] and |Rom03lj ) allow us to derive static arbitrage price bounds on a set of prod- 
ucts linked by a semigroup structure. The resulting constraints can be formulated as suc- 
cessively tighter linear matrix inequalities, hence we can compute increasingly sharp bounds 
on the solution to problem as solutions of increasingly large semidefinite programs (lin- 
ear programs on the cone of positive semidefinite matrices). Semidefinite programming has 
been the object of intensive research since the seminal work of |NN94j and several numerical 
packages (see for example SEDUMI by [Stu99] ) are now available to solve these problems 
very efficiently. 

The core of our argument is to substitute to the classical duality between the cones 
of probability measures and positive portfolios, the conic duality between positive definite 
functions on one hand and sums of squares on the other. These last two cones have the 
advantage of being numerically tractable and lead to exploitable formulations of the static 
portfolio super/sub- replication problems. 

A lot of work has been focused on arbitrage bounds in a dynamic setting, see |EKQ91|, 
EKQ95|, ALP95 j , and KS98 . among others. Work on the unidimensional static problem 
dates back at least to B1.78) (see also |LL00j ). both using the positivity of butterfly spread 



prices to preclude arbitrage. [BP02j studied these bounds together with second order moment 
constraints and proved the NP-Hardness of the multivariate problem Finally, in a 

previous paper |dEG03j . we focused on the interpretations of problem (JTJ as an integral 
transform inversion problem or a linear semi-infinite program, i.e. a linear program with a 
finite number of linear constraints on an infinite dimensional variable, and used the related 
theories to compute closed-form solutions for some particular cases and a linear programming 
relaxation for the general case. 

The paper is organized as follows. In section two, we describe the static market structure 
and start with a brief introduction on harmonic analysis on semigroups. Based on these 
results, we then derive necessary and sufficient conditions for the absence of arbitrage in the 
static market, formulated as semidefinite programs. Finally, in section three, we describe 
the conic duality between positive definite functions and sums of squares and use it to show 
how a super/sub-replicating portfolio can be constructed from the solution to the programs 
of the preceding section. 
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2 Static arbitrage constraints 



2.1 Market structure 

We work in a one period framework and suppose that the market is composed of cash and 
n underlying assets X{ for i — 1, . . . , n with x G R™ . We suppose that the forward prices of 
the assets are known and given by Pi, for i — 1, . . . , n, hence Wi is the Euclidean basis and 
Ki — for i — 1, . . . , n. In addition to these basic products, there are m + 1 basket straddles 
on the assets x, with payoff given by \w-[ +i x — K n+i \, i = 1, . . . ,m. Because a straddle is 
obtained as the sum of a call and a put, we get the market price of straddles from those of 
basket calls and forward contracts by call-put parity. 

We will note these payoff functions e,, for i = 0, . . . , m+n, with ej(x) = Xi for i — 1, . . . , n 
and e( n+ j){x) = \wfx — for j = 0, . . . , m. In what follows, we will focus on the Abelian 
(commutative) semigroup (§, •) generated by the payoffs ej(x) for i — 0, . . . , m + n, the cash 
1§ and their products. 

In this one period setting, we will look for conditions that guarantee the absence of static 
arbitrages, i.e. arbitrage opportunities that only involve trading today and at maturity, 
assuming that there are no transaction costs. 

2.2 Harmonic analysis on semigroups 

We start by a brief introduction on harmonic analysis on semigroups, for a complete treat- 
ment see |BCR84j and the references therein. Unless otherwise specified, all measures are 
supposed to be positive. 

Definition 1 A function p : § — >• R is called a semicharacter iff it satisfies p(st) = p(s)p(t) 
for all s, t G S and = 1. 

In |BCR84j an involution operation is defined on the semigroup (§, ■), here and in the 
rest of the paper we suppose that involution to be the identity, which means in particular 
that we take all semicharacters to be real valued. The dual semigroup of S, i.e. the set of 
semi characters on S is called §*. In this context, we call a function / :§^Ra moment 
function on § iff /(Is) = 1 and / can be represented as: 

f(s) = [ p(s)dv(p), for all s G S, (2) 

where v is a Radon measure on S*. 

When § is the semigroup defined in (j2.1J) as an enlargement of the semigroup of monomials 
on R n , its dual S* is the set of applications p x : § — » R such that p x (s) = s(x) for all s G § 
and all x G R n . The measure v is then assimilated to a probability measure on R n and the 
representation above becomes: 

f(s) = [s(x)] , for all s G S. (3) 

Our objective below is to find tractable conditions for a set of prices po, . . . ,p n+m to be 
represented as Fi u [\wfx — Ki\\ = pi for i = 0, . . . , n + m and some positive measure v. 
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2.3 The compact case 

In this section we assume the asset distribution has a compact support K. We treat the 
compact case independently as it is rather simple yet captures many of the key features of 
the general result. We begin by a few definitions along the lines of [BM84J and BCR84J. 
An absolute value on § is a function | • | : 8 — » R + satisfying 

|s 2 | < \s\ 2 , for all s G § 

and 

|ls| > 1- 

A function / : § — > R is said to be bounded with respect to an absolute value | • | iff there 
exists some M > such that 

\f(s)\<M\s\, forallseS. 

Furthermore, / is called exponentially bounded iff / is bounded with respect to some absolute 
value. Remark that if the measure v in has its support contained in the compact K then 
the moment function f(s) = JL p(s)dv(p) is bounded with respect to the following absolute 
value: 



\s\ K = supp(s) 
pgr 



for s G §. 



Definition 2 A function f : § — > R is called positive semidefinite iff for all finite families 
{s^ of elements o/§, the matrix with coefficients f(siSj) is positive semidefinite. 

We remark that moment functions are necessarily positive semidefinite. Necessary and 
sufficient conditions for the existence of a measure v in Q were derived in |HS90j . they 
were however numerically intractable. Here, based on the results in |BCR84j . |PV99j and 
|Rom03j . we look for exploitable conditions for representation (JHJ) to hold. 

Let a be an absolute value, the central result in |BCR841 Th. 2.6] states that the set of 
a-bounded positive semidefinite functions / : § — > R such that /(Is) = 1 is a Bauer simplex 
whose extreme points are given by the set of a-bounded semicharacters. Hence a function 
/ is positive semidefinite and exponentially bounded if and only if it can be represented as 
f(s) = L» pdu(p) with the support of v included in some compact subset of S*. 

Based on these results, we derive below a set of tractable necessary and sufficient condi- 
tions allowing a function / to be represented as in (j3J). For s,u in S, we note E s the shift 
operator such that for / : § — ► R, we have E s (f(u)) = f(su) and we let £ be the commutative 
algebra generated by the shift operators on §. Finally, we let (3 = sup x6i ^{^"^" m ei(x)}. 

Theorem 3 Suppose the asset distribution has compact support K and S is the payoff semi- 
group defined in \2.1)) . with [3 is defined as above. A function f(s) : S — > R can be represented 
as 

f(s) = E v [s(x)\, for all s E S, (4) 
for some measure v on K , and satisfies the price constraints in (QP if and only if: 
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(i) f is positive semidefinite, 

(ii) E e .f is positive semidefinite for i — 0, . . . , n + m, 
(Hi) ((31 — X]™iT-^ e i) / ^ s positive semidefinite, 

(iv) f(ei) = pi for i = 1, . . . ,n + m. 

Furthermore, for each function f satisfying conditions (i) to (iv), the measure v in repre- 
sentation is unique. 

Proof. The family of shift operators r = {{E ei }i = o,,.,,n+m, ifil — X/i=o™ -^i) } C £ is such 
that J — T G span + r for each T G r and span r = S, hence r is linearly admissible in the 
sense of [BM84, Corollary 2.5] or [Mas77J, which states that (ii) and (hi) are equivalent to 
/ being r-positive. Then, jMas77| Th. 2.1] means that / is r-positive if and only if there 
is a measure v such that f(s) = JL p(s)dv(p), whose support is a compact subset of the 
r-positive semi characters. This means in particular that for a semi character p x G supp(^) 
we must have p x {e,i) > 0, for i — 1, . . . , n hence x > 0. The set of r-positive semicharacters is 
then included in the nonnegative orthant and includes both the simplex {x > : ||x||i < /3} 
and K, hence / being r-positive is equivalent to / admitting a representation of the form 
f(s) = E u [s(x)], for all s G § with v having a compact support K C R". ■ 



2.4 The unbounded case 

The conditions derived in the last part do not describe all possible arbitrage free prices as 
they cannot account for unbounded asset distributions. Here, we use results from |PV99 
and |Rom03j to derive intrinsic characterizations of viable multivariate straddle prices. 

We note A(E>) the R-algebra generated by the functions Xs '■ §>* — ► R such that Xs(p) — 
p(s) for all s G §. By construction, Xs(p) — E s p(l§), and for a polynomial p G A(S) with 
p = J2k QkXg k an d f° r p G §* we have pp(s) = J2k Qkp(sgk) for all s G §. When § is the 
payoff semigroup defined in (j2.1|) . we naturally have Xs(Px) = s(x), for all x G R n , s G S 
and p G §*. 

We now note .4e(S) the R-algebra generated by A(S>) and 6 where 

0{p)= (l+> ]xM I , for all P G§*, (5) 



i=0 



we also note A(S, y) the algebra generated by A(S) and R[j/]. We first simplify the equality 
constraints on 2 n variables in |Rom03| Th. A] to recover an additive formulation as in 



PV99j. We begin by proving the following lemma. 



Lemma 4 The kernel of the algebra homomorphism $: 

A{S,y) -> A e (S) 



(6) 

p(p,y) ^ $p = p(p,6(p)) 



is the ideal generated by a £ A(S,y) such that o~(p,y) = y(l + X]i=L™ +1 Xe 2 {p)) ~ 1- 
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Proof. We adapt the proof of [PV99 , lemma 2.3] and let p G A(E>, y) be such that p(p, 9{p)) = 
0, we write p(p,y) = J2k1k(p)y k with q k G A(S). We have: 

p{p, y) = p(p, y) - p{p, 9{p)) = J2 k >o Qh(p) (y k - (d(p)) h ) 
= (y-(9( P )))l(p,y,(9(p)), 

where I is a polynomial. Let k = max{k : q k 7^ 0} and 

(m+n+l \ K 

1+ £ Xe? (p)J , for all peS*, 

we then have 

t (p)p(p, y) = ^(p, jOKp> y), (7) 

with r(p,y) G *4(S, y). The case k = is trivial hence we can assume k ^ 0. Using the 
fact that the polynomials r(z) and cr(z) have no common zeroes in C m+n+2 [z], Hilbert's 
Nullstellensatz (see [BCR98J for example) states that there must be f, a G C m+n+2 [;z] such 
that 

rf + 00 = 1. 

Multiplying this last identity by p yields, together with ((7|): 

p = a{rf + pa) 

hence the desired result. ■ 

The next proposition is adapted from the dimensional extension method in (PV99j Th. 
2.5] and [Rom03, Th. 4], to replace the exponential number of equality constraints in |Rom 03 . 
Th. A] with an additive formulation as in jPV99j . The function 9(p) is defined as in (J3J) and 
Ae(S>) is the R-algebra generated by -4(S) and 9. 

Proposition 5 With § being the payoff semigroup defined in ( POP ; let A be a positive 
semidefinite linear form on Ae(§) such that A (a^r 2 ) > for all r G and i — 1, . . . , n, 

then A has a unique representing measure v with support in R" and Ag(S) is dense in L 2 {y). 

Proof. We recall that the linear form A is positive semidefinite iff A(r 2 ) >, for all r G ^4e(§). 
As in |Rom03j . we define a bilinear form on r G Aq (S) by: 

{r 1 ,r 2 ) := A(nr 2 ), for all r 1} r 2 G A g (S) 

We let M be the set {r G ^4e(§) : A(r 2 ) = 0}. The bilinear form above then defines a scalar 
product on Ae(S)/N", and we note H the completion of this space. We define in TC the 
operators: 

Ti(r + J\f)= Xei r + Af, for all r G A (S)/Af and % = 0, . . . , n + m, 
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which are symmetric and densely defined in TC. We also define the operator (D(B),B) by: 



m+n 

2 



D(B) = M§)/M &ndB = J2 T r 



i=0 



The operator B is positive as a sum of squares of operators and, by construction, the domain 
D(B) is dense in H and invariant by B. Let r = Y^hLq 1 Xe 2 (p) an d r G Ae(S)/J\f, then u = rO 
is such that (1 + t)u = r, hence the operator / + B is bijective on D(B). This means that B 
satisfies the hypothesis of |PV99t Lemma 2.2] and is essentially self-adjoint. |Rom03| Prop. 
1] then implies that the operators T; for i = 0, . . . , n + m are essentially normal and that 
their canonical closures commute, meaning that there exists a common spectral measure H 
for the operators Tj for i — 0, . . . , n + m. With T = (Ti)i =Qt ... )n+m and r 6 Ae(§)/Af, we 
define the operator r(T) by: 

w+J\f i — ^ r(T)(w +Af) =rw +J\f. ^ ' 

With 7(x) = (X]r=o mx ?) > there is an element g of R 7 [x], the R-algebra generated by R[x] 
and 7(x) such that r(p) = q ((XeJ?.=o,...,n+m(p), for all p G §*. We then have: 

A(r) = (r(T)l, 1) = (g(T)l, 1) = / q(x)dH 1+NA+M (x), 

The homomorphism /: 

p(x) H-> /(p) =p((x ei (p))i=0,...,n+m,^(p)) 

satisfies the hypothesis of [Rom03, Lemma 2] hence there is a (positive) Radon measure v 
on such that: 

A(r) = / r(p)du(p), 



which, if § is defined as in section (|2.1|) . is also: 

A(r) = / r(x)dv(x) 



J HP 

Uniqueness and density follow from the argument in |Rom03j . Now, because the operators Tj 
for i — 1, . . . , n are essentially self-adjoint with A (x«r 2 ) > for r G ^4e(S) and i — 1, . . . , n, 
we know that the Tj are positive for all i. The spectral measure of Tj is given by Fi(X) = 
H (T i ~ 1 (X)) for all Borel sets X C R and Fi must be concentrated in R + for all i = 1, . . . , n 
hence the spectral measure H of T is concentrated in R™ and so is the representing measure 
v. ■ 

We can now formulate a general moment theorem that describes all the price systems 
that admit a representation as in ©• 
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Theorem 6 Let § be defined as in \2.1)) . A sequence f(s) : E> — > R is a moment sequence 
and can be represented as in 

f{s) = E v [s(x)\, for all s E S, 

for some measure v with support in R™, if and only if there is a sequence p(s, k) : (S, N) — > R 
such that: 

(i) p(s,0) = f(s) for all s E S, 

(ii) p{s,k) is positive semidefinite on (S, N), 

(Hi) p(eiS, k) is positive semidefinite on (§, N) for % = 1, . . . , n, 

(iv) p(s, k) = p(s, k + 1) - ESTX^s, k + 1) for all (s, k) E (S, N). 

Furthermore, the representing measure for sequence f is unique if and only if the sequence 
p is unique. 

Proof. First we show that conditions (i)-(iv) are necessary. With S, the payoff semigroup 
defined in (|2.ip . we recall that S* can be identified with R™, hence Xs{Px) = s(x), for all 
x E R+, s E S and p E S*. Suppose that / can be represented as: 



we let 



/ s(x)dv(x), for all s E S, 



m+n 

p(s,k) 



„ / m+n \ 

/ s(x) 1 + Y] e?(x) ) rfz/(x), for all (s, jfe) G (S, N), 
\ i=0 / 



which satisfies (i) and (iv) by construction, p(s, k) is then a moment sequence on the product 
semigroup ((§,•) x (N, +)) and as such must be positive semidefinite, hence condition (ii). 



Then, because for i = l,..,nwe have 



„ / m+n \ k 

I s(x) ( 1 + e?(g) ] ei(x)dz/(x), for all (s, fc) G (§, N), 
V i=o / 



m+n 

p(e»s, fc) 



we know that p(eiS, k) is a moment sequence for the measure ei(x)du, hence condition (iii). 

Conversely, let's assume that we are given a sequence p(s,k) satisfying (i)-(iv). We let 
*A#(S) and A(S,y) be the R-algebras described at the beginning of the section. We define a 
linear function A on A(B>, y) by: 



V j,k / j,k 
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and as in lemma we can define the following algebra homomorphism $: 



A(S,y) 

p(p,y) 



(10) 



whose kernel M has been computed in lemma EJ and ^(S) is isomorphic to the quotient 
A(fs>,y)/N. Condition (iv) implies that L{J\f) = and we can then define a linear form A 



with r G -4.6>(S>) and q G A(S,y). Because of (i)-(iv), the form A satisfies the hypothesis of 
proposition Eland has a unique representing measure v. m 



3 Price bounds and static hedging 

In this section, we show how the duality between the existence of a pricing measure and that 
of a replicating portfolio transposes into the moment framework described in the previous 
section. In particular, we detail how an optimal static super/sub-replicating portfolio can be 
constructed using the solution to the dual of to the moment problem in (fTJ). In particular, 
in a result that is consistent with the dynamic framework (see |ALP95j ). the replicating 
portfolio only involves options in the data set and no other option is needed to "complete 
the grid". 

3.1 Price bounds via semidefinite programming 

Here, we show how one can compute bounds on the solution of problem using a subset 
of the moment conditions imposed by theorem El These conditions cast (131) as a semidefinite 
program (see |NN94j or |VB96j ). which can then be solved efficiently using solvers such as 
SEDUMI by |Stu99] . 

3.1.1 Asset distributions with compact support 

As before, we note A(S) the R-algebra generated by the functions Xs '■ §* —* R such that 
Xs(p) — p( s ) f° r all s G § and p G 8*. For a polynomial p G A(S>) with p = YliQiXgi where 
gi G S, and for s G § we set 



With S the payoff semigroup defined in 1)2.1)1 . we recall that §* can be identified with R n , 
hence Xs(Px) = s(x), for all x G R n , s G § and p x G §*. This means that p G A(S) can be 
rewritten 



on A 9 (S) by: 



A(r) = L(q), where r(p) = q(p, 9(p)), for all p G §*, 
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We now recall the construction of moment matrices as in |CF00j and |Las01j . We adopt the 
following multiindex notation for monomials in .4(8): 

e a (x) :=e«»Wer(x)-e:r;W, 

and we let 

De = (1; CO; • • • > Cm+ri) e e l; • • • j e 0' • • • J e m+n) (H) 

be the vector of all monomials in A(§>), up to degree d, listed in graded lexicographic order. 
We note s(d) the size of the vector y e . Let y G R s ( 2d ) be the vector of moments (indexed as 
in y e ) of some probability measure v with support in R™, we note M d (y) G R_ s ( c( ) xs ( d ) ) the 
symmetric matrix: 

M d {y) i;j = (y e ) i (x){y e ) j (x)du{x), for i, j = 1, s(d) 

In the rest of the paper, we will always implicitly assume that y\ — 1. With /3(i) the exponent 
of the monomial (y x ) i and conversely, i(/3) the index of the monomial in y e . We notice 
that for a given moment vector y G TV^ ordered as in (|14j). the first row and columns of 
the matrix Md(y) are then equal to y. The rest of the matrix is then constructed according 
to: 

M d (y)i,j = Vi(a+p) if M d (y) itl = y i{a) and M d (y) ltj = y m . 
Similarly, let g G «4.(§), we derive the moment matrix for the measure g{x)dv on R™ 
(called the localizing matrix in [CFOOJ), noted M d (gy) G S^, from the matrix of moments 
M d (y) by: 

M d {gy) itj = / {y e ) i (x) (y e ) j {x)g{x)dv{x) 
Jul 

for i,j = 1, s(d). The coefficients of the matrix M m (gy) are then given by: 

9ayi(/3(i)+/3(j)+a) 

(12) 

a 

We can then form a semidefinite program to compute a lower bound on the optimal 
solution to Q using a subset of the moment constraints in theorem^ taking only monomials 
and moments in y up to a certain degree. 

Corollary 7 Let N be a positive integer and y G R^ 2 ^, a lower bound on the optimal value 
of: 

minimize po := E u [eo(x)} 

subject to 'E u [e i (x)] = Pi, i = 1, . . . , n + m, 

can be computed as the solution of the following semidefinite program: 

minimize y2 

subject to M N {y) >z 

M N (e jy )hO, forj = l,...,n, (13) 

M N ((P-Et™e k )y))hO 

U(j+2) = Pj, for j — 1, . . . , n + m and s G S 
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where s is such that i(s) < s(2N). The optimal value of (Q3J) converges to the optimal value 
of the original program as N — > oo. 

3.1.2 Unbounded distributions 

Here we work on the product semigroup (§, •) x (N, +). Its dual is the set of functions 
p x : (§, N) — ► R such that p x ((s, k)) = s(x)x k for all s E §, k E N and x E R n . As before, 
we note A(B>, N) the R-algebra generated by the functions Xs '■ (S, N)* — > R such that 
X(s,k)(p) — p{{ s ,k)) f° r all s e S and p E (§, N)*. With S, the payoff semigroup defined in 
(12. 1J) . here (S, N)* can again be identified with R n , hence X{s,k){.Px) = s(x)x k , for all x E R n , 
(s,k) E (S, N) and p x G (§, N)*. By construction, we have 

(X(s,k)) 2 = X(s*,2k), fo r all (s, k) E (S, N), 

and for a polynomial p E A(E>, N) with p = J2iQ.iXgi xki where (gi,ki) E (S, N), and for 
(s, Z) G (§, N) we set 

P((s> 0)0*0 = 5^gis(a:)^(x)x fcl+z , 

i 

for all x G R™. We adopt here the multiindex notation for monomials in A(S, N): 
e a := (e ,0r( ei ,0r---(e m+n ,0r- + "(l,ir- + " +1 . 

We then let 

y e = (1, (e , 0), . . . , (e m+n , 0), (1, 1), (e , 0) 2 , (e , 0)( ei , 0), . . . , (e , 0) d , . . . , (1, l) d ) (14) 

be the vector of all monomials in A(§>, N), up to degree d, listed in graded lexicographic 
order. We note s(d) the size of the vector y e . The matrices Md{y) and Md(gy) are defined 
as in the compact case above. 

We can again form a semidefinite program, this time using a subset of the moment 
constraints in theorem |3 taking only moments up to a certain degree. 

Corollary 8 Let N be a positive integer and y E H s( - 2N \ a lower bound on the optimal value 
of: 

minimize po := E u [eo(x)} 

subject to 'E u [e i (x)] = Pi, i = 1, . . . , n + m, 

can be computed as the solution of the following semidefinite program: 

minimize y2 

subject to Mn(v) b 

M N {{ ej , Q)y) y 0, for j = l,...,n, (15) 

En+m 
i=0 Vi(e?s,k+l) 

y(j +1 ) =pj, for j = l,...,n + m and (s, k) E (S, N) 

where (s, k) are taken such that i(s, k) < s(2N). The optimal value of hlty converges to the 
optimal value of the original program as N oo. 
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3.2 Static hedging portfolios and sums of squares 

We let here S C A(S) be the set of polynomials that are sums of squares of polynomials in 
.4.(8) , and V the set of positive semidefinite sequences on S. The central argument of this 
paper is to replace the conic duality between probability measures and positive portfolios: 

p(x) > <^ J p[x)dv > 0, for all measures v, 

by the conic duality between positive semidefmite sequences V and sums of squares polyno- 
mials S: 

(f,p) > for ah> G £ iff / G V. 

for p G A{E>) with p = J^iliXst and / : § — > R having defined (f,p) = J2iQif( s d- Tlie 
previous section used positive semidefmite sequences to characterize viable price sets, in this 
section, we use sums of squares polynomials to characterize super/sub-replicating portfolios. 
From the initial price problem written in terms of straddles: 

minimize p := J Rn e (x)dv(x) 

subject to J Rn ej(x)<iz/(x) = pj, i — 1, . . . , n + m, (16) 
/ Rn dv(x) = 1, 

in the variable u, a positive measure on . We can form the Lagrangian: 

(x) - Xn+m+l 1 dl>(x) 



n+m, „ / n+m 

L(u, A) = Xn+m+i + ^2 XiPi + ( e Q (x) - ^ Xiei 

i=l >/ R™ V i=l 



with variables v and A G R n+m+1 . We obtain the classic dual as a portfolio replication 
problem: 

maximize \ n +m+\ + **Pi n 7 \ 

subject to e (x) — YliZ™ ^i e i( x ) ~ X n+m+ i > 0, for all x G R™ . 



in the variable A G R n+m+1 . 

Unfortunately, the problem formulations above are numerically intractable except in 
certain particular cases (see BP02J and [dEG03 ). On the other hand, as we have seen in 
the previous section, the conditions of theorem0]turn problem into an infinite dimensional 
semidefinite program which can be relaxed to produce tractable bounds on the solution of 
(JTJ). Here, we detail the accompanying duality theory to exhibit a static hedging portfolios 
corresponding to these bounds. 

We can assume without loss of generality that the payoff functions {ej(x)}j = o,..., m +n, 
together with the cash 1§, are linearly independent. Then [BCR84, Proposition 6.1.8 and 
Theorem 6.1.10] hold and we can form a dual to the cone of positive semidefinite functions 
on S as follows. For p G «4(§) with p = £\ qiX Si and / : S — > R with: 

(f,p) =5^?»/(si), 
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BC R841 Theorem 6.1.10] states that E is the polar cone of V for the above bilinear form, 
in other words: 

(f,p) > for allp G £ iff / G V. 

We can use this conic duality to compute a dual to program (JT3J). Considering the compact 
case for simplicity, Corollary (JJJ) states that the initial pricing problem: 

minimize p := Ej,[e (:r)] 

subject to Ej,[ei(a;)] — Pi, i — 1, . . . , n + m, 

is equivalent to the following (infinite) semidefinite program: 

minimize yi 
subject to M(y) y 

M(ejy) >z0, for j — 0, . . . , n + m 

Vti+2) = Pi> for j = 1, . . . , n + m 
l/i = l- 

in the variable y : S — ► R. We can form the Lagrangian: 

L(y, A, q) := y% + (1 - yi)\ n +m+i + Eji™ (Pj ~ V(j+2))Xj - (y, ?o) 



ELT ( e ^> <?;> - ((P ~ EL™ e k)y, Qn+i) 



or again: 



L(y, K q) ■= 2/2 + (i - yi)X n +m+i + E"=r (Pj - 2/y+2))Aj - (y, <? ) 
- EJ= m e i?i) - (0 - EST efc)g n +i) 

in the variables y : § -> R, A G R n+m+1 and ^- G £ for j = 0, . . . , (n + 1). We then get the 
dual portfolio problem: 

maximize YJjtT Pj\ + ^n+m+i 

subject to e (x) - E"=i™ \ e j( x ) - K+m+i (18) 
= Qo(x) + J2]=T <lj( x ) e j( x ) + e k (x))q n+1 (x) 

in the variables A G R n+m+1 and ^ G £ for j = 0, . . . , (n + 1). 

The key difference between this portfolio problem and the one in (|17j) is that the (in- 
tractable) positivity constraint e (x) — ES™ \ e i( x ) ~ A n+m+ i > in (jT7j) is replaced by 
the tractable condition that this portfolio be written as a combination of sums of squares 
of polynomials in ^4(§). Such combinations can be constructed directly from the dual solu- 
tion to the semidefinite program in (|13p. hence a numerical solution to the program in (JT3~j) 
provides both a price bound and an accompanying portfolio. 
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